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IntroducCton; The need for an Intiiractlve Biosphere 
In General Circulation Models. 

Importance of Evapotranaplratlon . 

The sensitivity of weather and climate to land-surface evapotransplra- 
tlon is difficult to determine from observations, but has been revealed in 
experiments with numerical general circulation models • (Sec Shukla and 
Mintz, 1982; and the paper by Mlntz, 1982, in which eleven such experiments 
are reviewed.) In each experiment two calculations were made, which started 
from the same initial atmospheric state and had the same sea surface tempera- 
tures and sea-ice extent, but with land-surface conditions that produced 
different evapotransplratlons . In every case the result was a large differ- 
ence in the calculated precipitation, temperature and motion field of the 
atmosphere. 

A decrease in latent heat transfer from the land to the atmosphere 
is approximately balanced by an increase in the sensible heat transfer. 

The sensible heat transfer, however, warms the air within the relatively 
shallow planetary boundary layer, and it does so locally in space and 
time. The latent transfer, by contrast — if realized as heating through 
the convective condensation process — warms the air in the free atmosphere 
up to the tropopause level. Moreover, because of water vapor advectlon, 
the realization of the heating may be at nome distant place and a later time. 
It is the difference in the vertical distribution of the two kinds of 
heating, as well as the possible shift in their horizontal distributions, 
which makes the general circulation (the thermally-forced large scale 
atmospheric circulation) sensitive to land-surface evapotransplration. 
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Importance of Vegetation . 


Because numerical weather predictions and predictions of climate 
anomalies with general circulation models will be sensitive to the land- 
surface evapotranspiration, great care must be used in the way in which 
the evapotransplration is calculated. 

In almost all existing general circulation models (see the review 
by Carson, 1982), an ’’open bucket" formulation is used to calculate the 
transfer of water vapor from the land to the atmosphere. The level of the 
water in the bucket is lowered when evaporation is larger than precipi- 
tation; and the level is raised when precipitation is larger, up to the 
point at which the bucket overflows and produces "runoff". Over some 
broad range in the level of the water In the bucket (which varies somewhat 
from model to model), the rate of evaporation is taken as equal to (or 
nearly equal to) the evaporation from a free water surface; and only when 
the water level is low is the evaporation rate made less than that from a 
free water surface. But this is hardly the way in which water vapor is 
transferred from the land to the atmosphere in the real world. 

A more realistic representation of how water is exchanged between 
the land and .atmosphere is shown in Fig. 1 (from Rutter, 1975). Here 
the atmosphere is insulated from the water in the soil by a vegetation 
layer; and the largest part of the water transfer to the atmosphere is the 
transpiration of water which the roots of the plants take up from the. 
soil. Another large part of the precipitation (and under some vegetation 
and precipitation conditions the largest part) never enters the soil at 
all, but is Intercepted by and stored on the surface of the plants and 
from there evaporated into the air. The third transfer process, the 













6 


original page 
OF POOR QUALITY- 


evaporation of water from the pores of the solli is generally much smaller 
than the other two transfers. 

The v^ater storage capacities on the surface of the plants, within the 
plants and within the plant root zone of the soil, are of the order of 1, 

10 and 100 ram, respectively. Inasmuch as both the Interception loss and 
the transpiration are of the order of a few mm per day, the characteristic 
recycling time for the intercepted water is of the order of a few hours, 
and for the soil moisture store it is of the order of a month. 

The 100 mm of water which the vegetation can remove from the soil is 
only a small fraction of the total water in the ground. Similarly, of the 
total of about 30 mm of water vapor held in an atmospheric column, only 
about 3 mm , rut of about the 10 mm in the planetary boundary layer, can be 
removed from the atmosphere by the convective precipitation process. 
(Reducing the relative humidity in the PEL by more than about 30% shuts 
off the convection.) Thus, if not renewed, water vapor in the atmosphere 
can provide only about 1 day's worth of space-averaged convective precipi- 
tation. Soil moisture, however, even when not renewed, can maintain the 
average rate of evapotransplration for about a month. It Is this which 
makes the soil moisture store more important than the atmospheric water 
vapor store for numerical weather predictions and predictions of monthly 
and seasonal anomalies of climate. 

Vegetation exerts physiological and morphological controls over evapo- 
transpiration. Because transpiration is associated with photosynthesis, 
it is usually confined to the daylight hours; unlike the open bucket formu- 
lation, which can produce appreciable evaporation at night. Moreover, 
when water is available, the rate of the bucket evaporation is limited only 
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by Che aCmospheric factors of radiation, wind speed and air humidity j and 
Che daytime rate, especially In semi-arld regions, can be very large. But 
with plants there is also a daytime limit to the transpiration rate, because 
there is a constraint on the rate at which water can flow through the plants. 

The morphology (the vertical structure) of the vegetation has a large 
influence on the rate of interception loss. With tall (forest) vegeta- 
tion, the elevated and dispersed surfaces of the leaves are more strongly 
ventilated than are the surfaces of the leaves in short (herbaclous) vegeta- 
tion. When the leaf surfaces are wet, therefore, the interception loss 
from Call vegetation is much greater than that from short vegetation. But 
when the leaf surfaces are dry, the greater ventilation of tall vegetation 
is generally compensated by its larger storaatal resistance; so that (under 
the same atmospheric conditions) the transpiration losses for tall and 
short vegetation are not very diffetent. 

The interception loss rate from short vegetation is not much larger 
than its transpiration rate; but with Call vegetation, the interception 
loss rate can be 5 to 10 times larger than the transpiration rate. Even 
though the vegetation surface storage capacity is small, Che time-integrated 
interception loss from tall vegetation can be very large. 

Fig. 2 is an example of Che mean annual water and energy balances of 
short and tall vegetation, as derived from measurements in two adjacent 
catchments which have nearly the same atmospheric conditions. The total 
evapotransplration loss from the forest catchment is more than twice as 
large as that from the grass-covered catchment* With the grass-cover, 58% 
of the net radiatlonal heating of the surface is used for evapotransplration, 
and 42% for the sensible heating of the atmospheric boundary layer. But 


with the foreat cover, the energy used for evapotranspiration exceede the 
radiatlonal heating by 15 %, and this results in a removal of sensible 
heat from the boundary layer (there is u negative Bowen ratio, --0.13). In 
the forest case it was possible to measure the two components of the evapo- 
transpiration separately, and the interception loss was 1.7 times larger 
than the transpiration. Although the atmospheric conditions are about the 
same for the two catchments, the net radiatlonal heating, R^, is about 
lOX larger with the forest-cover than with the grass-cover. This is not 
only because the forest is darker and absorbs more of the Incident solar 
radiation, but also because, as a consequence of the larger evapotrans- 
piration rate, the surface temperature of the forest is lower and emits 
less infrared radiation. It is clear, from numerous examples of this 
kind, that vegetation exerts a large control over the water and energy 
exchange between the land-surface and the atmosphere. 
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Fig. 2 Measured mean-annual water and energy balances of adjacent 
grass-covered and forest-covered catchments, in central Wales, U.K. (from 
Calder and Newson, 1979; Shuttleworth and Caldcr 1979). 
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At Goals of the Proposed Research. 

The objective of the proposed research Is to construct a model biosphere 
which will produce realistic simulations cf the water and energy transfers 
at the earth's surface. This will be done by constructing a numerical model 
of the earth's vegetation cover which embodies the principal morphological 
and physiological factors that control the water and energy transfers. 

Morphology 

The morphological factors (those which depend on the form and structure 
of the plants) will be the leaf area density (the total area of the leaf 
surfaces per unit volume of space) as a function of height. Lq[Z]; and 
the root length density (the total length of live roots per unit volume of 
space) as a function of depth, RTq(Z). 

In the initial version of the model biosphere, Lp[Z] and RTd[Z] will 
be prescribed for each grid area of the GCM as a function of the time of 
the year, as obtained from phenologlcal observations. 

In the second stage, Lp[Z] and RT[)[Zl will be made interactive with 
the model calculated atmospheric conditions (and soil moisture), so that 
aperiodic drought and extremes of heat and cold will also affect these two 
parameters on the phenologlcal time scale. 

In the final stage, the different vegetation formations themselves 
(l.e., rainforest, seasonal forest, woodland, desert; grassland) and not 
just their phenologlcal changes will be made interactive with the atmospheric 
conditions (and soil moisture.) 

PRECEDING PAGE CLANK NOT FILIICO 
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PhyaioloBv 

The principal physiological factor will be the stomatal resistance^ 
rg(. (the impcdonce to water vapor transfer from the saturated cavities 
within the loaves to the air outside*) rgj. will depend not only upon 
the atmospheric variables, but also upon other laternal lasistancoe which 
affect the water flow from the soil to the leaves via the roots and stems* 



B» Overview of the Proposed Dloaphcre 


As described In Section B»2» the vegetation of the model biosphere 
will have two or more discrete canopy layers and two root layers. The two 
important vertically continuous vegetation properties that will be dis- 
cretized are the leaf area density as a function of height, and 

RTd[Z 1, the root length density as a function of depth. 

B . 1 . 1 Governing equations for transpiration and water uptake by roots . 

For the purpose of illustration, we show here the governing equations 
for the simple case where there is one canopy layer and one root layer. 

For the definitions of the symbols, see Section C, p 39. 
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If the vegetation is represented as a single, continuous, transpiring 
surface of unit area per unit ground area, the governing equations for the 
energy transfers can be written: 


- Cfl) pCp 
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Given 6a, Tg, Ug, Sg and RLa from the output of the atmospheric part of 
the general circulation model; £ from the heat budget of the soil (see 
section Iv) ; and from the solution for the water uptake from the soil, 
as described below, we can solve these six equations for the six unknowns: 
Ef, e^, T^, H, tg, rgt* 


IS 


(1) As shown scheraatlcalXy in Fig, 3, the main part of the transpira- 

tion loss is water which evaporates from the walls of the raesophyll cells 
that surround the sub-stomatal cavities in the leaves of the plant and then 
diffuses to the atmosphere through the stomatal openings. [The rate of 
water vapor transfer across the cuticle surface of the leaves » or across 
other parts of the plant surfaces, is usually one to two orders of magnitude 
smaller.] The rate of the water vapor transfer, from its origin in the 
sub-stomatal cavities to a given reference level in t!je atmosphere, is 
given by Eq. 1, where is the vapor pressure in the sub-stomatal cavity, 
and is the vapor pressure at the reference level in the atmospheric 
boundary layer, rgt is the bulk resistance to the diffusion of water 
vapor through the stomatal opening, and r^ is the resistance to the 
diffusion of water vapor fJrom the vegetation surface to the reference 
level. With eg given, Eq. (1) has 4 unknowns: Ex» Ojl, rgt» Cg. 

(2) can be taken as the saturation vapor pressure, e*, as given 
by the Clausius-Clapyron equation for the temperature of the leaf, Tjj. 

T^ is assumed to be constant throughout the leaf (whose thickness is 
typically about 1 mm) and represents an additional unknown. 

(3) T;j is obtained from the diffusion equation for the sensible 
heat transfer from the surface of the leaf, at temperature to the 
reference level in the atmosphere, at temperature, Tg, where it is 
assumed that r^ is the same for the sensible heat diffusion as for the 
water vapor diffusion; which adds the unknown, H. 

(4) H, the sensible heat transfer, is obtained from the total vertical 
energy transfer, where the sura of the three terms within the parentheses 

on the right, in Eq. 4, is the net radiation flux, Rjj. Sg and Rlq 
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are the dovmward components of the solar and longwave radiation, as given 
by the atmospheric part of the GCM, a Is the albedo of the vegetation 
surface, and oTn** Is the upward longwave radiation, which depends on the 
temperature of that surface. G is the sensible heat transfer through the 
vegetation and into the soil, as given by the heat budget for the soil 
(section iv) . 

(5) The aerodynamic resistance, r^, depends on the vertical distri- 
bution of the leaf area density, on the wind speed at the atmospheric 
reference level, and on the thermal stability of the air, (Tjj -Tq) , (see, 
for example, Goudrlaan, 1977). Typically, with wind speeds in the atmos- 
pheric boundary layer of a few meters per second and mid-day unstable lapse 
rates, r^^ is of the order of 0.5 s cm”^ for herbacious vegetation (0.1 to 

1 ra high) and of the order of 0.05 s cm”^ for forests (10 to 20 ra high). 

(6) rgi;, the bulk stomatal resistance reaches the minimum value for 
the plant species, rgj-^o, when the stomates are fully open. Typically, 
when the leaf area index is about four, r_,. . is of the order of 0.5 s cm'"^ 
for herbacious plants, and somewhat larger, of the order of 1 .0 s cm~^ , 
for trees. 

The stomates close when the solar radiation, S^, drops below a 
critical value they close when the relative humidity of the air, 

RHa » RHtea, Ta), falls below a critical value, RHa,crlt» 
close under extremes in the leaf temperature, Tj^. In addition, the 
stomates close when the water potential in the guard cells (ij)j^) that 
surround the stomates, falls to a limiting value, 'j'Jl,crit» (What 
controls the magnitude of i)i£ is discussed below.) When the stomates 
are closed, the transpiration is only through the cuticle of the leaf. 
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Fig. 3 Schematic representation of transpiration through the 
and cuticle of a single leaf (after Hillel, 1971). 


Btomates 
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whose resistance to the water vapor transfer is typically about two orders 
of magnitude large., than r^j., or about 50 to 100 s cra”^ when equals 

four. 

All these responses are .interpreted as the evolutionary result of the 
plant’s effort to maximize photosynthetlc production (which requires the 
simultaneous presence of shortwave radiation, water, carbon dioxide, and 
sclutes in the leaf chloroplast cells), to maintain a solute transport sys- 
tem, to avoid heat death, and to conserve water in the soil moisture store 
for future use. The first three demands are In direct conflict with the 
last, the conservation of water, which gives rise to the need for physio- 
logically regulating the flow of water through the plant. 


quality 

Water Uptake by the Roots; 
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The conservation of mass requires that the water vapor which moves 
away from the cell walls of the sub-storantal cavities bo replaced by a 
flow of liquid water toward the walls on the other side. 

If, for the purpose of Illustration, we neglect the small time rate of 
change of the water stored in the plant, we have 


Ecp - Fp “ Ur , 

where E«p is the rate of transpiration, Fp is the rate of liquid water flow 
through the plant, and Ur is the rate of water uptake by the roots. 

When there is no d’.vergene.e in the water flow through the plant (no 
change in the plant water storage), we can write the governing equations 
for the rate of water flow as 


C'i's “ 'I'i) 


Fp “ 

(7) 

g(rB + Tc + r^) 


ts - 'I'stll 

(8) 


rs(0, RTdIZ)] 

(9) 


(10) 

Fp^crit) • 

(11) 


Given the boundary condition, Fp « E-p, where E-p is obtained from the solu- 
tion of equations (1, 2.... 6) above; and given 0 from the solution for 
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the soil water budget (section iv), we can solve these five equations for 
the five unknowns; (Jig, rgj r^. and r^* 

There is an upper limit to Fp, which is believed to be the rate at 
which the water movement through the fine channels in the xylem (the water 
conducting elements of the stem) changes from laminar to turbulent flow. 

This limit, Fp^crit* typically of the order of 1.5 mm/hr (when expressed 
as the area averaged flow). Whe.n Fp^^Pij^j. is reached, there is an abrupt 
and large Increase in the xylem resistance, r^^*, and, as Eq. (7) shows, 
for a given soil water potential, \|Jb, this will produce an abrupt and large 
fall In the leaf water potential, Itself has a critical value (of 

the order of -25 bars) which is when which the water pressure in the stomatal 
guard cells can no longer maintain their turgor. When the cells collapse 
they constrict the stomatal opening, and E<j ■ Fp is then limited to the 
value Fp^erlf 

Fig. 4 is an electrical analog of the water transfer pathway through 
the soil-plant-atraosphere system (when Ej “ Fp “ Uj^) . The figure shows 
how the flow is controlled by the fixed resistance, r^, and the four 
variable resistances, rg, r^, rgg and ra. 
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Fig. 4 Electrical analog represej>CaCion of the potentials, and re- 
sistances to flow, of water vapor and liquid water In the soll-plant- 
atraospher<; system. 
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ii. Rainfall Interception and Interception Loss . 

For the single layer vegetation model, the governing equations for 
the Interception loss and Interception gain, when the canopy Is completely 
wet, can be written: 






[fia - eal pCp 

' T - T • ■ ■ 

Ta Y 

{ 0, when Wj ■ 0 

1^8 “ ®al pCp 

ra Y 


f (Wi) , When Wi > 0 


when Wj > Wx ^ax 


( 12 ) 


Wi max Is the maximum amount of water that can be held on all the 
upper leaf surfaces. 


dWi 

P (1 - p) - Ex - D [Wxl . (13) 

dt 


The Interception loss and the Interception gain are highly dependent 
on the vegetation morphology for two reasons: 

1) The aerodynamic Interaction between plant and atmosphere; As 
Indicated in section B.1.1, tall, 'rough,' surfaces, like pine forests, 
maintain a relatively turbulent aerodynamic regime In and around their 
upper crowns. Within this well ventilated volume, the transfers of vapor 
and sensible heat are typically an order of magnitude faster than at an 
equivalent height above 'smoother' surfaces, like a grass cover, under the 



same meteorological conditions, (r^ ~ 0,05 s cm~^ as against 0.5 s cm~^) 

11) The geometry and extent of Intercepting surfaces; Vegetated 
regions commonly present 3 to 8 ra^ of surface area to the atmosphere for 
every square meter of ground area « 3-8). Ohviously, the larger 
the surface area the greater the amount of rainfall that may be held on 
the surface for later re-evaporation, and the smaller the amount that 
drains off to reach the soil moisture store. Typically, the Interception 
capacity of vegetation is of the order of 1 to 3 mm depth of water. The 
geometrical arrangement of the canopy is important too - some coniferous 
trees maintain a large Intercepting surface well above the level of the 
theoretical momentum sink. 

The rate of interception loss depends very greacly on the rainfall 
regime. For a given time-averaged rainfall which is made up of high in- 
tensity, short duration episodes, there will be a comparatively low rote 
of interception loss. This Is because the interception capacity is quickly 
reached after the onset of storms, allowing the remaining rainfall to be 
transmitted to the soil surface. Conversely, if the same time-averaged 
rainfall is made up of long-duration, low intensity storms, the interception 
loss will be larger, reaching Its maximum if the Interception capacity is 
not reached (l.e. if the water does not accumulate on the leaf surfaces to 
a level sufficient to Initiate significant leaf drainage). 

As Indicated in section B.2, our intention is to model the water and 
heat transfer processes for the several observed vegetation formation types 
In a realistic fashion. This, however, will produce correct values of 
the interception loss only if the space-time variation of the rainfall is 
correctly represented. General circulation models produce half-hourly to 
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hourly amounte of rainfall as output; but as this Is the average over a 
large grid area (of the order of (300 km)^), If used In an unadjusted 
form it will almost always produce an excessively high and unrealistic 
interception loss. It will be necessary, therefore, to relate the GCH 
produced grid area average rainfall to a realistic space variation of the 
rainfall on the sub-grid scale, using parameters extracted from the GCM 
simulation (which indicates, for example, whether the rainfall is of con- 
vective or large-scale upgllde origin) and also climatology. Satellite 
based observations of regional surface temperature and soil moisture varia- 
tions (Atlas and Thiele, 1981) may be able to provide information about 
typical tlrae-area-intenslty distributions in the different regions and 
seasons, from which we would construct the area-intensity functions for 
the GCM rainfall calculation. A possible solution liss in representing 
the spatial variation of the amount of water held on the vegetation surface 
by wave functions. A similar wave representation of the rainfall intensity 
would be superimposed to obtain the spatial variation of the rate of change 
of the water held on the surface. Another possibility would be to assume 
that the rainfall at a representative point within the grid area will be 
equal to the grid-area averaged rainfall when both are averaged over a 
time interval of the order of, say, 12 to 24 hours; and then to let the 
half-hourly values of the rainfall intensity at that point vary within the 
chosen time Interval as a function of the GCM parameters and climatology. 



ill. Snowfall Intarceptioii and Dtsposal 


The partition of snow into surf ce evaporation and melt water, which 
contributes towards soil moisture recharge and runoff, is of great Importance 
for the surface energy balance and hydrology of a large part of the Northern 
Hemisphere continents. The energy balance at the snow surface can be 
modelled using derivations from the equation set outlined in section B.1.1. 
Major differences between the fluxes of latent and sensible heat with and 
without snow, under the same atmospheric conditions, arise from variations 
in the surface albedo and hence in the net radiation, and from the value 
of the aerodynamic resistance, r^* [A review paper by Male and Granger, 

1981, discusses these and other facets of snow surface processes in detail.] 
The energy available for evaporation and/or snowmelt in forested regions 
is characteristically much greater than that over bare or grass covered 
areas (Leonard and Eschner, 1968) due to the smaller albedo of the exposed 
parts of the forest canopy. This produces an appreciable transfer of 
heat, in radiative and sensible form, from the exposed parts of the tree 
canopy to the snow covered surfaces. This effect is etjhanced in areas 
^d.th evergreen vegetation. The use of a multilayer vegetation model 
coupled with a layered soil model ia the most realistic way to simulate 
the energy exchange processes involving snOw' cover. 


Iv. Water and Heat PudRet o£ the Soli 
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For almuXatlng the water and heat budgets of the 0 oil» we shall follow 
one of the more physically realistic treatments (see, for example, the reports 
of Gurney and Camlllo, 1982, Camlllo end Schmugge, 1982) t In these treatments, 
the fluxes of heat and moisture are coupled as both depend on the gradients 
of temperature and moisture. This requires the simultaneous solution of 
the equations: 


% “ D0/ii\ - Vir\-Kg 
\dz J \dz } 



where q 0 and qj^ are, respectively, the fluxes of moisture and heat; 
lh’, 0 ,vap dlffuslvlty of heat, water, or water vapor^* and 

Is the thermal or hydraulic conductivity of soil. It should be noted 
that Dq and have both liquid and vapor contributions; 


Dq " ^0,liq + D0 vap 

(15) 

Dt - Dx,iiq + Dt 

,vap 

Equation (14) describes the dependence of the heat fluxes on the vertical 
gradients of 0 and T. The time dependence of 0 and T is given by the 
continuity equations: 


dt 

dt 

where Cg 


Vdz / 

- 

^dz j 

Is the specific heat of the soil. 


( 16 ? 
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Various numerical techniques may be oppllcd to solve the equation set 
over finite time steps. It Is usual practice to prescribe a deep soil 
temperature and moisture content as lower boundary conditions and an energy 
balance model at the soil surface to calculate *10, 0 and T ai the alr- 

ooll Interface. The soil Is then divided Into a number of strata, which 
are taken as Internally homogenous with uniform values of ond 

The values of T and 0 for each stratum are used to calculate the depth 
dependent derivatives at each time step. 

The soil moisture store will be represented by three zones: a small 

near-surfacc store, from which not only root uptake but also direct evapora- 
tion can take place; an Intermediate bulk soil moisture store, which Is 
mainly drawn on by the root uptake; and a deep store, from which only 
capillary rise can bring water toward the surface. The last terra can be 
significant on the seasonal time scale. 

The hydraulic and thermal properties of the soli will be taken from 
the global data set prepared for the GLAS GCM by Lin and Alfano (Alfano, 
1981). 

It is Impossible to model overland flow and soli interflow, In a 
direct way, without resort to an exceedingly fine grid size. Instead, we 
will use the ‘lumped* catchment analytical model of Beven and Klrkby (1976). 
This makes use of functions which relate Che total soil moisture storage 
to the size of the catchment contributing area, (which Is the saturated 
area bordering stream channels), thus allowing a direct calculation of the 
different components of runoff: overland flow. Interflow and base flow. 

It may be possible to use the geographical data relating to grid square 
soil type, topography and stream density to parameterize the necessary 


functions 
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B»2 The Multilayer Vegetation Model, 

1 ) Definition and Advantages of a Multilevel Model ■ 

Definition; 

In general, a multilayer model can be thought of as an elaborate 
extension of the short set of equations given in Section B.1.1, The vege- 
tation canopy Is represented by a vertical series of discrete plates, 
each of which exchanges sensible and latent heat with the Immediately 
surrounding air and through which the fluxes must pass on the way to the 
free atmosphere (see Fig. 5). At the same time, the temperature of each 
plate depends on the net radiation It absorbs (which Is a function of the 
temperatures of all the plates and the soil surface), and on the temperature 
and vapor pressure of the air in the surrounding (canopy) air space. It 
is apparent from Figure 5 that, given; (1) the air temperature (T^) 
and vapor pressure (e^) above the canopy, and their equivalents below 
the soil surface; (2) the values of the intermediate resistances and (3) a 
means of distributing the shortwave and downward longwave radiation among 
the plates, the equation set representing the energy budgets of all plates 
will reduce to n expressions with n unknowns, (where n is the number of 
leaf layer plates and the unknowns are the leaf temperatures.) A similar 
model will represent the layers of the soil. The combined system will 
represent the steady state energy balance of the vegetation and soil, and 
may be further extended to Include non-steady state processes - such as 
the drying out of Intercepted water which involve a dependence on the 
time variation of canopy and soil temperatures and water contents. 
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Fig. 5 Schematic representation of sensible and latent heat transfer 
using a multilayer model. The flow of water from soil to leaf layers has 
been omitted for clarity. 
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Advantages; 


The multilayer model Is more realistic and Is easier to validate than 

a unllayer model* The principal points to note are: 

a) Vegetation morphology: A multilayer model can be constructed In such a 

way that the leaf area density, LptZj, and the leaf area Index, 

are* directly represented by the size and spacing of the leaf layer 
'plates’ . 

b) Energy exchange processes: The exchanges of energy (In radiative, 

sensible and latent form) between the different layers of canopy, and 
between those layers and the soil, can be modelled. It need not be 
assumed that the multiple latent and sensible heat sources correspond 
to the momentum slnh, as with a unllayer model. 

c) Parameter correspondence with nature: The parameters of the model can 

be made to correspond directly to the physical and physiological proper- 
ties of each vegetation type. 

d) Interception loss and snowmelt; These processes cannot be modeled 
realistically using a unllayer treatment. Sellers (1981) demonstrated 
that a single-layer model consistently underestimates the rate of 
Interception loss. 

e) Validation; A multilayer model not only generates the latent and 
sensible heat fluxes, but also profiles of leaf temperature, air 
temperature, vapor pressure, soil temperature, and soil moisture 
potential. By comparison, the unllayer model outputs a single value of 
’surface’ temperature - a rather nebulous quantity when considering 
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that a tropical ralnforeat has been observed to have a mld~a£ternoon 

temperature difference of as much as 7° C between the top of the canopy 

and the soil surface (Pinker, 1980; and personal communication, 1982). 

11. Distribution of the Observed Vegetation Forma t ions. 

de Laubenfels (1975) has made the following classification of the 
earth’s vegetation formations, whose distributions are shown In Fig, 6; 

Region with tall trees whose crowns form a continuous canopy, and 
below which there Is a continuous understory of shorter trees: designated 

Rainforest . 

Region with trees whose crowns form a continuous canopy, but below 
which there Is a discontinuous understory: designated Seasonal Forest . 

Region with trees whose crowns do not form a continuous canopy, but 
where the total vegetation cover Is continuous: designated Woodland . 

Region where there Is a discontinuous cover of plant growth, so that 
large areas of bare ground are exposed: designated Desert . 

Region of continuous ground cover of herbaclous plants (of which 
grass is the predominant form): designated Grassland . 

Rainforest, Seasonal Forest, Woodland, and Desert are "undisturbed" 
vegetation formations. Grassland Is the existing formation type where 
seasonal forest and woodland have been disturbed by fire and by grazing 
(which Inhibits recovery after fire.) 





Seasonal forest 

4 I -I i / J 



Fig, 6 Observed vegetation formations (from de Laubenfels, 1975) 
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In the Initial version of the proposed biosphere, the prescribed 
distribution of the vegetation formations will not change with time* 

Eventually, we would like to have the vegetation formations change in 
response to forcing by the atmospheric part of the GCM* But Inasmuch as 
the time scale of natural succession Is of the order of 10^ to 10^ years 
(see, for example, Loucks, et« al. , 1981) this will be an essential require** 
ment for the biosphere only when a fully Interactive ocean 1s made part of 
the general circulation model. 



iil) Specification of the Multilayer Model 
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There arc two distinct requirements for the specification of the 
multilayer model. Firstly; the morphology of the vegetation community 
must be realistically represented by the size and spacing of the discretized 
leaf plates (see Fig. 5)» and by the values of the boundary layer resistances, 
canopy air space rcslstonces and root cortex resistances. The above canopy 
aerodynamic resistance, r^, must also be properly related to the gross 
morphology of the vegetation community. Secondly; the physiology of the 
plant must be correctly reflected in the functional model of stomatal 
response and its interaction with the leaf water potential, xylem resistance, 
stem flow, and soil moisture potential. 

Taking the morphology of the vegetation community first; Figs. 7A, 

Q, C illustrate the process whereby measurable characteristics of the 
vegetation are transformed to model parameters. Fig. 7A shows the structure 
of a tropical rain forest in its natural state. Fig. 7B shows the corres~ 
ponding leaf area density and root length density as functions of height. 

In the first version of the model, the seasonal variations in these vegeta- 
tion parameters will be prescribed as functions of the time of the year. 

Later, the phenolgical changes in LpCZ) and RTd(Z) will be made dependent 
upon the model derived atmospheric variables and soil moisture. 

To obtain the areas of the two plates that will represent the vegeta- 
tion (see figure 7C), the leaf area density is Integrated with height; 

Lo - / Ld(Z) • dZ (17) 

where is the leaf area index (m^/ra^), and Zg and Z^ are, the heights of 
the bottom and top of the storey. 


Abstract 


The purpose of adding an interactive bio^iphere to the GLAS general 
circulation model is to make the calculation of land-'surface evapotransplr- 
ation and sensible heat flux more realistic and, thr;rcfore, more accurate. 
This is important because sensitivity experim^mts vlth general circulation 
models have shovjn that the large scale fields of rainfall, temperature and 
motion of the atmosphere arc highly sensitive to the transfers of latent 
and sensible heat at the land surface. 

Water and energy tronsfers at the land surface depend on vegetation 
morphology and physiology. In the proposed model biosphere, the vegetation 
morphology will control these transfers through the leaf area density os a 
function of height, Ld(Z), and through the root length density as a 
function of depth, RTjj(Z), as expressed in a discretized form for each 
of the model grid areas. Lj)(Z) will influence (a) the aerodynamic resistance 
to the latent and sensible heat transfers; (b) the transmission, absorption 
and emission of radiation energy by the canopy and underlying ground; and 
(c) the interception, evaporation and throughfall of rain and snow. The 
vegetation physiology will control the transfers of water and energy through 
the stomatal resistance, the xylem resistance of the stems, and the root 
cortex resistance. The stomatal resistance will be a function of the leaf 
water potential, short-wave radiation intensity, leaf temperature, and 
humidity of the air. 

In the initial formulation, Lp(Z) and RTp(Z) will be prescribed as 
functions of latitude and longitude and the season of the year, as known 
from ecological observations in the various vegetation formations. Later, 
the phonological changes of Lp(Z) and RTp(Z) in the deciduous forests 
and grasslands will be made interactive with the atmospheric variables 
and the soil moisture ns calculated by the GCM. Finally, the vegetation 
formations themselves will be made interactive with the atmospheric and 
soil moisture conditions. 

The realism of the model biosphere will be evaluated (1) through 
short period one-dimensional comparisons with measured values of the local 
atmospheric forcing and measured values of the vertical fluxes of water 
and energy; (2) through long period comparisons of simulated and observed 
catchment water budgets; and (3) through comparisons of the simulated and 
observed surface temperatures, surface albedos, snow cover, live biomass, 
and the water balances of the large river basin'^, when the model biosphere 
is forced with atmospheric variables taken from the IGGE level 3-B and 
2-C data sets. 
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The height of Che plate la taken as the center of gravity of the plate, 
Z„, where: 

Zj^ Zf 

/ Ld(Z) • dZ - / Ld(Z) • dZ 

Zb z„ 

Following the work of Goudriaan (1977), the bulk aerodynamic parameters 
of roughness length (Zq) and zero plane displacement (d) may be determined 
from a numerical analysis of the reduction in shear stress with depth in 
Che canopy; parameters which depend on Che leaf area density and drag 
coefficient of leaf elements. The parameters Zq and d arc used to determine 
the value of Tq, the above canopy resistance, by 



(19) 


where k is von Karmart's constant, h is the height of the crop and is the 
reference height (This equation holds for neutral conditions only; but Thom 
and Oliver (1977) describe a version which accounts for the effect of non- 
neutrality). A further analysis by Goudriaan (1977) provides the extinction 
of wind speed as a function of leaf area density. 

The conductivity of mass, heat and momentum in the canopy air space 
is given as; 






( 20 ) 


where is the local mixing length of the canopy air space (a function of 
the size and spacing of shoot elements' and is the relative turbulence 
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Fig. 7: A. 


B. 


C. 


Tropical rainforest vegetation; continuous canopy of tall 
trees above continuous under-storey. 


af area density * seen area density * 

ot length density HrHIQ , as functions of height for the 
onical rainforest shown In Figure A. 


Electrical resistance analogue of the tropical rainforest 
represented In Figure B. 
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;lntcnalt;y » a faccor char Increases from 0<3 ac Che Cop of plane canopies Co 
0*8 aC Che base, and appears Co be InvaclanC wich wind speed* 

r^ , Che canopy air space rcsiscancc between one leaf layer and Che nexC 

Sg 

is given by; 

H-l 

re " / * dZ (21) 

The leaf laminar boundary layer resistance, tj, , will be taken as a 

i 

function of leaf area index, local wind speed, and a shape/shcltcr factor, 

Xs5 

rb ” * u. * Xs> (22) 

i (z) (z) 

Several researchers, notably Goudriaan (1977) and Allen and Lemon (1976), 
have proposed different machematical analyses to determine Xs a function 
of leaf size, shape and orientation. 

The radiation absorption and transfer characteristics of the vegetation 
are determined by the canopy structure. A simple model of radiative transfer 
was proposed by Ross and Nilson (1971): 

Sa,, , - Sa,„,U(w) + exp(-f(3) • Lj,)] (23) 

lOJ 


where Sq Is the shortwave radiation flux emerging from under a cumulative 
(l*(y) 

leaf area index radiation intensity above the canopy, 

f(w) is a scattering function, and f(8) la an extinction coefficient that 


varies wich the solar angle, 3* 
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The root water uptake models of Cowan (1965), Hillel (1977) and others 
describe the movement of water from soil to root cortex as a radial flow 
towards Individual root elements across a gradient of moisture potential - 
the difference between i|)g, the soil water potential, and the root 
water potential* The plant parameter of Importance Is RTd[ZJ, the root length 
density (mm'’^), as a function of depth. Cowan (1965) originally proposed 
the expression; 

D /r^\ (24) 

3t 7 3r \ 9r / 


where D Is the value of s'oll moisture dlffuslvlty, and r Is the distance 
from the root. Given that the roots have a typical radius of r^ and that a 
root extracts water from a cylindrical volume of soil of radius r 2 » where 
r 2 ■ (1 / ir.RTolZl ] equation (24) can be solved with the boundary condi- 
tions of 


/ 


fits - 

6r 




d0 - r^2] 

— at r ■ tj 

dt 2r^ 


0 


at r •* T 2 


where dWg/dt is the rate of reduction of soil water concentration. 

The transfer of water from root cortex to leaf mesophyll involves 
aspects of plant morphology and physiology. The xylem elements, which 
conduct water up the stem are serially linked, elongated plant cells with 
narrow apertures at their end-to-end junctions. Poiseuille analyses of the 
flow In such elements (see, for example, Denmead, 1976) predicts a critical 
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flow velocity above which the flow regime changes from laminar to turbulent. 
This greatly increases the resistance imposed by the xylcm and has the 
effect of lowering 

The end of the liquid transfer pathway is the leaf raesophyll coll. The 
leaf water content, , is determined by: 

dWji 1 

«— P— ^ ^ 

dt g 

Leaf water potential, a critical variable in the calculation of 
stomatal resistance, Is a function of Wj^, 

^ - £(v»t) . 

The particular form of stomatal response is dependent upon plant 
speclec, where each species can be thought of as being near perfectly 
adapted to survive in its natural environment. In the main, stomatal response 
is a function of shortwave radiation intensity, leaf water potential, 
vapor pressure, temperature and leaf age. Closure of the guard cells is 
brought about by a reduction in their turgor; this may be Induced by a slow 
hydropassive process, whereby water is abstracted from the guard cells by 
strictly pViyslcal processes - such as evaporation or osmotic flow - or by a 
fast hydroactive process. It is believed that the latter, which can cause 
complete stomatal closure within seconds, is controlled by the transport 
and metabolism of ions into and out of the guard cells under the control 
of a plant hormone, Absclsslc acid (ABA). The local concentration of ABA 
is determined by (see figure 8). 

It is not proposed to model the biochemical and detailed biophysical 
responses of the leaf, although this has been done by some researchers 


I*! - tr ~ 
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Fig.e 


Dependence of Tq); 

1975) 


on \|»jj and CO2 concentration (From Raschke, 
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(aec, for example the report of Penning do Vries; 1971); but rather to 
make the storaatal resistance a function of the leaf water potential and 
atmospheric forcing. Thus; 

’’^8t(®a) 

^Bt • == (25) 

f('I'l) • f(ea) • f(^) 

Examples of the Individual functions in equation (25) are given below. 


Light 


rst(Sa) - Ai/(A 2 ~ A 3 • (26) 

Denmead and Millar (1976) 
where Aj^ is a species dependent constant. 

Leaf water potential 


f (t|ijj) Ml- exp(-A 4 • 6 \(/;^) 

" tpji - >l'i»crit (27) 

0 < f(i{<;i) < 1 Turner (1974). 

where tjl,crit *^he value of below which the stomata close completely. 
Vapor pressure 


\e^) " ~ CeMTal - ^1 


0 < 


f 

(ea) 


< 1 


( 28 ) 


Jarvis (1976) 
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Temperature 


- ’’bX'^u - 

Aj - 1/(T„ - T^)(T„ - 
A7 - (Tu - To) /(To - Tb) 


0 < f(Ta) < 1 


Jarvis (1976) 


(29) 


where, Tb and Ty are the temperatures above which and below which the 
Gtomata open, and To is the temperature at which £( Tg ) “ !• 
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B«3 ComparlsonB of the Model Performance with Observations* 

The checking of the model's performance xd.ll be done In three phases, 
each of v;hlch corresponds to a further level of horizontal Integration. 

These are briefly stated belovf, starting at the lowest level of Integration. 

I) Field Measurements. 

The complete vegetation-soil model for each vegetation type xdll be 
driven offline using local mlcrometeorologlpal data. There are a number of 
case studies reported in the literature xjhere a short time series of simul- 
taneous observations xvere made of heat and vapor fluxes, temperature and 
humidity profiles, and mlcrometeorologlcal conditions above the canopy. 

Use of these data sets should ensure that the characteristics of the various 
vegetation types are transformed to the model in a physically reasonable 
xvay. 

II) Water Balances in Catchments. 

Here, the comparisons are made over a time period of a year or more 
for an area where long records of meteorological and streamflow observations 
exist. This will allow a checking of the simulated water balance components 
against the observations for the different vegetation types. It should 
be noted that in both (1) and (11), no allowance Is made for horizontal 
differences In meteorological conditions: the tests are essentially one- 


dlmenslonal 
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ill) FGGE Observatlona for the Globe» 

Uolng a 4~dlmcn8lonal dynamic assimilation of ground based and satellite 
based observations of the atmosphere during the FGGE year (1979) > (^hc GLAS 
GCM Is producing a time series of the atmospheric state variables, Tq, 

C|j and Uj,; and of the radiation transfers and R^a, averoged 
for the model grid areas and at half-hourly time lnt>ervals. In addition, 
the GIJIS model generates grid-area averaged rainfall, [P]; which will be 
converted Into representative hourly local rainfall Intensities, as Indi- 
cated In section B.l.li. The observed dally precipitation at 30,000 stations 
arc available on tape for the FGGE year, and can also be used for Input 
( 10,000 of these are In the U.S. and have hourly records.) 

This information will be used to force the model biosphere and, thereby, 
obtain a time series of the global fields of Ei, Et, H, and oTs**} Wj and Wg; 
and of the drainage, Dgj for the FGGE year. 

The model derived OTg*^ , and its spatial and temporal variation as a 
function of Wx, snow cover, etc., (where Tg is the radiation temperture 
of the earth’s surface) can be compared with satellite measurements of Tg. 
This will provide a quantitative evaluation. 

Ex can be compared with the observed distribution of actively growing 
vegetation. C. J. Tucker has shown that growing and dying herbacious 
vegetation can be determined from Landsat observations. This, however, can 
only be a qualitative evaluation. 

Ex, H and Wx will be difficult to evaluate. Wg^l, the moisture in the 
upper root zone, whose depth Is of the order of 10 cm, may perhaps be 
compared with satellite microwave measurements of the moisture in the 
uppermost 5-15 cm of the soil in the desert and grassland regions (Schmugge, 
et.al., 1980 ). 
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Dg can be compared with the mcaaured river flows during 1979, which 
are now being documented, on tape, as part of the FGGE Level 2~B data act. 

•• It la the above (phenologlcally Interactive) version of the biosphere 
that should be used when making weather .predictions and predictions of monthly 
and seasonal climate anomalies with the general circulation model# 
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B,4 Numerical Simulation of the Vegetation Formations* 

In the biosphere described above , the phenological changes o£ the 
vegetation formations will Is interactive with the atmosphere and soil 
moisture, but the vegetation formations themselves will be prescribed (as in 
Figs • 6 A”E) ■ 

If the above goal of biosphere modeling is successfully accomplished by 
the beginning or early part of the third year, we shall undertake the next 
step: which is to model and simulate interactive vegetation formations. 

Here we will use empirical expressions which relate the climax vegetation 
type to the atmospheric and soil moisture conditions and thereby (starting 
from some initial state of the earth’s vegetation cover) derive the distribu- 
tion of the "undisturbed" vegetation formations: the forests, woodlands and 

deserts. Replacement by grassland would occur wherever atmospheric and 
soil moisture conditions make the forest and woodland vegetation tinJe’-dry; 
the assumption being that some mechanism for setting the dry vegetation on 
fire (lightning, dewdrops acting as burning lenses, or incendiary man) is 
always present. 

A substantial extension of the above procedure would be to make use of 
the ecological processes which govern the growth, aging and succession of 
vegetation types. Starting from a given Initial state, we would calculate 
the long term changes in the vegetation formation parameters. 

Validation 

Validation can be made, on the simplest non-trlvial level, by forcing 
the biosphere with the non-interactive FGGE atmospheric data set. Here, we 
would assume that this one year data set represents the climate that produced 
the observed vegetation formations. 
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A higher level of validation would be to Inltallze the vegetation 
formations with the FGCE data set and then let the model atmosphere and 
biosphere operate In the interactive mode (with some large compression of the 
successlonal time scale.) Xn doing thlS| the presently observed ocean 
surface temperatures would be prescribed* Where the simulated and observed 
vegetation formations agreed, the model would be taken as being correct: 
where they disagree, we would conclude that either the model was not 
correctly formulated o£ that the presently observed vegetation formation 
Is not In Its natural equilibrium state. 

Applications. 

Atmosphere-biosphere fully Interactive; ocean surface temperatures 
prescribed. 

Besides the second level of validation Indicated above, there will be 
relatively few us‘es for a fully interactive atmosphere-biosphere GCM when 
the ocean surface temperatures are prescribed (or when the ocean temperature 
Is prognostic only In Its thin boundary layer.) With given ocean temperatures, 
we are limited to the CLIMAP type of calculation. For example, we can 
prescribe the paleontologically derived global ocean surface temperatures 
of 18,000 years ago, and derive both the atmospheric state and the vegetation 
formations of 18,000 years ago. 

Ocean~atmosphere-blosphere fully Interactive. 

Eventually wa shall have a global ocean model which is fully interactive 
with the atmosphere and which correctly simulates the heat transport and heat 
storage in the deep ocean. (Fpr a review of the near-current state of the 
art of ocean modelling, see Mlntz, 1979.) When that is In hand, there 
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be raany Interesting and important climate simulation and sensitivity experiments 
that can be made with the fully Interactive ocean-atmosphcre-blosphere GCM; 
some* concerning natural variations in climate ^ others concerning anthropomorphic 
Influences on climate. 

An example, in the first category, would be the simulation of climate 
when the orbital parameters of the earth with respect to the sun are changed 
(l.e., a numerical test of the Mllankovitch hypothesis on the Ice ages.) 

In the second category, which concerns man's Influence on climate, per- 
haps the most Important use of the fully interactive ocean-atmosphcre-blosphere 
model would be to study the changes In the entire climate system (Including 
the changes in the stores of carbon In the ocean, atmosphere and biosphere) 
as a result of the burning of fossil hydrocarbons, or of man’s raodlflcation 
of the living biomass (through desertification and deforestation, or conserva- 
tion and afforestation.) 


C. Notation 




c 


P 


c 


Q 


d 


DtWj.) 

ea 

e* 


ejl 

h 

Et 

f(ea) 

f(w) 

f(0) 

HH) 

V 

‘ p,crit 

g 

G 

h 

H 

I 

*ni 

k 


Bpccloa dependent constants 

specific heat of air (J kg"^ *0”^) 

(•• 1*01 at s.ttp*) 

specific heat of soil (J tn~^ '*0“^) 

zero plane displacement (m) 

dlffualvity of heat, water or vapor In the soil (m a 

drainage rate (Kg ra * s ) 

vapor pressure (mb) 

saturation vapor pressure (mb) 

vapor pressure In the substomatal cavity (mb) 

Interception loss rate (Kg ms*) 
transpiration rate (kg ms*) 
vapor pressure component of storaatal resistance 
air temperature component of stomatal resistance 
scattering coefficient of canopy 

radiation extinction coefficient as a function of solar angle 
leaf water potential component of stomatal resistance 

"■9 ••1 

flow of water through the plant stem (Kg ra""^ s * or mm hr ) 
limiting value of F (Kg s“^) 

r 

acceleration due to gravity (m s“2) 

heat flux to ground (W m“2) 

height of the crop (m) 

sensible heat flux to atmosphere (W m“^) 

relative turbulence intensity 

von Karraan's constant 
(» O.Al) 


PRECEDING PAGE BLANK NOT FILMED 
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^»v,h/ conductivity of raomentura, vupor and hoot in the air (ra a 


KO 

Vz) 

lm 


P 

P 

r 

>^a 

""h 

"c,o 

'"at 


(z) 

thermal conductivity of soil (W m“^ °C“^) 
hydraulic conductivity of soil (m a”^) 
leaf area denaity (ra^ m“^) 
local mixing length (m) 
leaf area index (m* m *) 
number of leaf layers 
throughfall coefficient 
rai .;all rate (kg m"2 a~l) 
radius (m) 

aerodynamic resistance (s m~^) 

leaf laminar boundary layer resistance (s ra“^) 

canopy air space resistance (a m“^) 

root cortex resistance (a) 

root cortex resistance per unit length of root (a m*'*') 
soil resistance (s) 
storaatal resistance (s ra"^) 


*^at(Sa) 


x,o 


R 


La 

% 

rTqCz) 

fa 
s 

o) 


component of storaatal resistance dependent on shortwave 
radiative Intensity only (s ) 

xylera resistance (s) 

xylem resistance per unit height (s m*"^) 

longwave radiation incident on the surface (W m“^) 

•-2 

net radiation (W m ) 

root length density (m m ) 

solar radiation incident on surface (W m”^) 

solar radiation incident above canopy (W m ) 

air temperature (°C) 
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T|, ccmperaturo below which BComaces close (**C) 

Tj deep soil temp era ture (**C) 

T)i loaf temperature (*C) 

Tg temperature at which otomatal roslatance factor, f(Tg) - 1, ("C) 

Tg surface temperature (”C) 

Ty temperature above which stomates close (*’C) 

u wind speed (ra s”l ) 

flow of water through the plant root system (kg ra'*^ s”'^ or mm hr”^) 
leaf water content (kg m"^ or mm) 

Wj water held on canopy surface (kg m“^ or mm) 

z height or depth (m) 

Z|j height of bottom of canopy (m) 

Z (1 reference height (ra) 

zi height of 1th layer above ground (m) 

Zm height of leaf 'plate' (m) 

Zg roughness length (m) 

Zf height of top of canopy (m) 

a albedo 

y psychrometrlc constant (mb “C"^) « 0.646 at s.t.p, 

X latent heat of vaporization (J kg"l) “ 2.501 at s.t.p. 

jjjjl leaf water potential (bars) 

'll!. root water potential (bars) 

ilfg soil water potential (bars) 

p density of air (kg m^^) 

(• 1.292 at s.t.p.) 

fl soil water content (kg m”^ or m^m''^) 

Xg leaf shape/shelter factor 
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